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ABSTRACT 


The Electrochemical Machining ( ECM) is one of the most 
widely used unconventional machining process. However, the 
potential of this highly capable process has not been fully 
exploited due to the inherent complexities of the process, 
like , simultaneous occurrence of two phase fluid dynamics, 
unsteady state heat transfer, mass transfer, thermodynamics, 
electrostatics and electrochemistry between moving boundaries* 
This complicates the problem of anode shape prediction in 
EQyi many folds. Majority of the models available for anode 
shape prediction are based on simplified assumptions. 

In the present work, a modified anode shape prediction 
model has been proposed which predicts the anode profile in 
stagnation, front, transition, side and stray current zones. 

This model accounts for the variation in different process 
parameters during machining. Experimental data has been used to 
ascertain the accuracy of the model. Comparison of analytically 
predicted and experimentally obtained anode profile reveals a 
reasonable agreement between the two. The effects of the 
assumptions made and the experimental errors and their possible 
contribution in causing the deviations between theory and 
experiments have been discussed. 
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Conventional ECM tool design is not accurate enough to 
account for different complex phenomenon taking place in the 
domain of interest. A method has been suggested for determi- 
nation of the shape of the tool (cathode) (specifically 
accounting for the effects of stray current zone) to produce 
the required work ( anode) shape, while machining under 
specified conditions , The tool design model has been further 
extended to predict the tool shape which will produce a spike 
of specified size in the stagnation zone. This model can be 
used for analysing the problems of tool design for external 
shaping by ECM, This method accounts for variation in different 
process parameters. The cathode shape prediction problem has 
been attempted for electrochemical drilling with a bare tool* 
Tools have been designed for producing profiles obtained 
experimentally* Tool design has also been carried out for 
straight sided holes with the top surface flat, as well as for 
a reasonably small spike .in the stagnation zone. Comparison 
of designed and experimental tool shapes reveal good agreement* 


t 
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CHAPTER I 

INTRODUCTION AND LITERATURE SURVEY 

1,1 ELECTROCHEMICAL MACHININO ( ECM) 

Advanced technological developments in industries like 
aerospace, nuclear etc, have been accompanied by development 
of materials which are exceedingly difficult to machine. These 
materials have been developed to meet the demand of high strengl 
and heat resistance. These rapid developments in the field of 
materials have given an impetus to the iranufacturing technology 
to develop, modify and discover newer technological processes 
with a view to achieve results that are far beyond the scope of 
existing conventional manufacturing processes. This has led to 
the development of new techniques of machining such as ECM, EDM 
USM, AJM, EBM etc. These unconventional processes are capable 
providing effective solutions to the problems imposed by the 
increasing demand for accurate machining of high-strength- 
temperature-resistant ( HSTR) alloys, the requirements of parts 
with complex geometries, and materials so hard as to defymachi 
ing by conventional methods. 

Of all the unconventional machining processes. 
Electrochemical machining is the one that is most widely used, 
ECM is a process of metal removal from an electrically conduct! 
workpiece by controlled and accelerated anodic dissolution. 

High velocity electrolyte flows between anode and cathode 
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(Figure 1,l) subjected to a small d.c, voltage and a small 
interelectrode gap (IBG) is maintained between them. The 
flowing electrolyte helps not only in allowing the high rate 
of metal dissolution, but also carries away the reaction 
products and heat ge’nerated during the electrochemical 
reaction. 

The advantages of ECM process lie in its ability to 
machine very hard metals without causing any tool wear, 
high metal removal rate, no structural damage to. the workpiece 
and capability to machine complicated shapes with good surface 
finish. The various areas of application [l ] of ECM are drilling, 
boring, turning, cavity sinking, milling, grinding etc. Although, 
the applications of ECM are many, its capabilities have not 
been fully exploited due to lack of clear understanding of 
complex nature of the process of metal removal and non-existence 
of an efficient tool design methodology,. 

The metal removal rate (MRR) in ECM depends upon a large 
number of interrelated parameters. Most of the machining 
parameters such as electric field distribution, electrolyte 
con djctivit y ( K) , temperature rise (At), valency of dissolution 
(Z) , electrolyte flow velocity (U), etc. are varying simultane- 

ously and continuously along the electrolyte flow path with 
machining time. Further, the interaction of these parameters is 
of complex nature which makes it further difficult to understand 
the exact mechanism of metal removal. Moreover, the electrolyte 
in the 103 is a mixture of liquid electrolyte, gases evolved 
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during electrochemical reaction and precipitates. In majority 
of cases, it is difficult to evaluate the flow parameters due 
to the effects of electrolyte starvation and cavitation. Thus, 
ECM is a fairly complicated phenomenon since it involves 
simultaneous occurrence of two-phase fluid dynamics, heat 
transfer, mass transfer, thermodynamics and electrochemistry 
betv/een moving boundaries. 

Further, the accuracy of machining is governed by how 
accurately the tools can be designed, vhich in turn depends upon 
the accuracy of anode shape prediction model. With this in 
view different models available in literature for cathode and 
anode shape prediction have been discussed below, 

1 .2 TOOL DESIGN FOR ECM 

The tooling problem is of great concern because it is 
usually linked with high cost, A direct method for deciding 
upon the tool shape is almost impossible for all but the 
simplest of work shapes. In practice, a 'trial and error’ 
method usually has to be relied upon, with repetitive machin- 
ing runs having to be made to obtain the correct tool geometry. 

The designing of tools for ECM can be classified under 
two heads. First is the determination of shape of tool 
together with optimum machining conditions necessary to 
produce the required work shape [2 1 , The second aspect 
is a practical one, which is concerned with rreking the tool 



of an appropriate material, fixing it in the machine, 
connecting it to the power supply and passing an adequate 
supply of electrolyte between tool and workpiece. The 
design requirements for some of these aspects might turn 
out to be conflicting and therefore may call for certain 
modifications in the tool geometry. Here, only part of 
the first aspect of tool design has been considered, i.e,, 
designing the tool shape for producing desired work profil< 
ECM deals with computation of tool shape which under specil 
machining conditions would produce a work having a prescril 
shape or profile. Inversely, the problem of anode shape 
prediction deals with determining the work shape obtainabli 
from a tool of known geometry when machining is performed 
under specified machining conditions [4,5], The workpiece 
profile is usually ccraputed from the equilibrium gap. The 
lEG is not uniform in ECM and it is influenced by a number 
of parameters such as electrolyte flow rate, viark material 
microstructure, change in valency of dissolution of the 
work-material during machining, type of electrolyte, currei 
density, stray current attack etc. The influence of all 
these continuously varying parameters is not well known, 
which makes the anode shape prediction difficult. 

In general, the practical methods like 'trial and 
error’ methods for obtaining the tool shape are expensive, 
time consuming and inaccurate. Several attempts have been 
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made to develop a comprehensive model of the mechanics of 
metal removal in ECM, Some of the models for anode shape 
prediction and cathode shape design are discussed below, 

1 .3 COS0 METHOD 

This method was first proposed by Tipton [6,7] for 
the computation of equilibrium gap for any given operating 
conditions. This method excludes the considerations of the 
mode of electrolyte flow, overpotential, variation in electrolyte 
conductivity, heat transferred to the electrodes etc. This 
method can be applied to the electrode regions where the 
electric field can be assumed to be nornpl to the electrode 
surfaces. In this method equilibrium workshape is computed 
corresponding to the tool whose profile has to be approximated 
by several planar sections inclined at different angles, say, 

G as shown in Figure ( 1 ,2) , 

For plane parallel electrodes with their surfaces 
normal to the direction of feed, the equilibrium gap, Y^, is 
given by 


K ( V ~ A V) 
^ f 


( 1.0 


where, K is the electrolyte conductivity, (V- V) is the 
effective voltage, F is the Faraday's constant, is the 

mass density of the work and f is the feed rate of the 


tool 


6 


The equilibrium gap between any section of cathode 
surface and the corresponding surface of the anode which 
is parallel to it is given by, Y /CosQo The angle '9' is 
measured between normal to the anode surface and the 
direction of the cathode feed. The cathode shape for a 
given anode surface can be computed from the following 
equations , 

Let Y = f C x) represent the anode surface. Then any 
point, sayA(x, y), on the anode surface has a corresponding 
point, say B(x^,y^) on the cathode surface. The gap width 
between A and B then becomes Y^/Cos9 (Figure 1,2), We have. 


Y-Y^ = AB CosQ = Yg 

a nd 


(1 ->3) 


X-j-X = AB SinQ = Y^ tan9, 
or 

• (I-'*) 

By substituting the values of X and Y from equations , 
(1.3) and (1,4) in equation (1,2) gives the equation for 
the cathode surface. 

The ’CosG’ method is generally applicable to the 
regions where the local radii of curvature of the anode and 
cathode surfaces are large compared with the equilibrium 
gap. This method excludes the effect of different process 


parameters on predicted anode profile. Workpieces having ' 
sharp corners and complex- shapes cannot be analysed with 
this method* 

1 ,4 NOMCGRAPHIC APPROACH 

K’dnig [4] and Ganesh fs] suggested the nomographic 
approach to evaluate the equilibrium anode shape, K<5nig [4] 
has prepared a nomogram for calculating the side gap for 
known values of equilibrium gap, tool radius and bare tool 
length. Ganesh [sjhas developed separate nomograms for 
evaluating the equilibrium gap and tool feed rate. Tool 
feed rate is calculated for the given current and area of 
cross-section of the tool. The equilibrium gap is then 
evaluated from the thus calculated feed rate and the given 
voltage. Nomograms facilitate planning of ECM operations, 
but since they have been prepared with many simplified 
assumptions they have not gained much popularity. Moreover, 
they are applicable for specific machining conditions and 
work-tool combinations,, 

1 .5 CONTINUITY METHOD 

Lawrence [9] developed a graphical technique for 
tool design by solving the field equation in the gap. 

But this procedure gives low accuracy and is tedious. He 
subsequently adopted numerical methods to solve the field 
equation using digital computer. In this method, starting from 
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a defined equipotential work shape boundary, a model of 
electric field is constructed, where each new equipotential 
surface is a possible tool shape which will produce the 
original work shape. This method rests on the assumption 
that Ohm's law applies across a small element of electrolyte 
and that the current flowing between a pair of flux lines in 
an electric field is constant. It does not take into account 
the variation in electrolyte properties in the IB3 along 
the flow path, 

1 ,6 CATHODE SHAPE PREDICTION BY DRAWING CURVILINEAR SQUARES 
IN THE IB3 

Tipton [10 ] has developed a method for determining 
the tool shape by drawing curvilinear squares between the 
work and tool boundaries. The work surface is defined numeri- 
cally in XY plane, Y is taken as feed direction and coordinates 
are spaced by constant X increment equal to a fraction of 
the equilibrium gap. The work boundary is considered as 
equipotential with each coordinate representing the end of 
a current flow line. In equilibrium, the normal current 
density being constant, a constant current increment 
flows between each of the adjacent coordinates on the work 
boundary since they are equally spaced in X-direction, By 
using each adjacent pair of coordinates, it is possible to 
draw the curvilin^r squares and to construct the first row ’ 
of potential field distribution between work and tool as shown 
in Figure 1,3. If the X coordinate spacing is then 
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applying the above process repetitively will result in 
th 

n equipotential which is also the tool surface. Each 
intermediate surface is also a possible tool surface for 
equilibrium gaps as Y /n, 2Y /n, 3Y /n, 4Y /n, etc, 

G € G © 

This method is error prone, especially iMien the 
process is repeated several times, since it is based on 
differentiating the input data for drawing curvilinear 
squares, and when repeated any random error in the input 
data gets magnified. Complete failure of calculation 
occurs when two current flow lines cross each other. In 
such a case, the previous equipotential could be used as a 
tool with the fractional equilibrium gap associated with 
it. This model does not account for variations in the 
process parameters, 

1 .7 NUMERICAL METHODS 

The emergence of fast digital computers has led to 
the application of numerical techniques to simulate the ECM 
process. In ECM, during machining, local current density 
causes the change in the shape of workpiece. Further, 
the current density in turn changes with the changing work- 
shape, To determine the current density, it is necessary 
to solve Laplace equation for the field distribution in 
the lEG. The three popular methods used to determine the 
potential distribution are discussed below^ 
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1.7,1 FINITE DIFFERENCE TECHNIQUE (FDT) 

Tipton [6 ] and, Hopenfeld and Cole [11] applied 
FDT to determine potential distribution in the IE3, 
assuming constant conductivity and temperature. 

Figure (1.4) shows the tool work surfaces drawn 
in square meshes containing a set of grid points of general 
coordinates ( i, j). The grid points on the cathode boundary 
are assigned zero potential values, and those on the anode 
boundary are set to some known values. The initial potential 
at the grid points within the IB3 region are set by linear 
interpolation along the vertical grid lines between anode and 
cathode. These values must be progressively adjusted until 
they satisfy a finite difference equation corresponding to 
Laplace's equation in the region between the boundaries. 

For a point 'O’ located in a mesh of spacing h(i,j), 
the finite difference equation corresponding to Laplace's 
equation becomes. 



Using equation (1,5), the potential at *0’ can be 
obtained as. 


= 




( 1 . 6 ) 
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In some cases, all points on the anode and cathode 
boundary ina y not lie on the grid points* In such cases, 
i t is customary to relax the potential by linear inter- 
polation, Nanayakkara and Larsson [12] have suggested 
the use of irregular grids along with regular grids to 
account for complex geometries, 

1.7.2 FINITE a^BlENT METHOD (FEM) 

In majority of cases, the lEG is complex shaped. 

As a result, the cathode and anode boundaries cannot be 
matched accurately using square meshes, which introduces 
further approximations. In order to overcome this problem, 

Jain and Pandey [13] suggested the application of FEM to 
analyse the complex shaped lEG, 

Jain [3] has developed FEM models, FET-11 and FET-22, 
for the prediction of anode profiles in case of plane parallel 
machining. The model FET-11 is based on use of Ohm's law 
for the calculation of current densities in the 103, and FEM 
has been applied for one dimensional problem to obtain 
temperature distribution v/ithin the 183. In the model FET-22, 
Laplace equation has been solved to obtain potential distribution 
w.Hhin the 183 from which current densities are calculated, 

Jain also extended model FET-11 (and named it as SGFET-11) for 
predicting the anode profile in case of rectangular and 
cylindrical deep hole drilling operations. All these models 
are capable to predict the anode shape accurately, by taking 
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into consideration, the variations in different process 
parameters, Murugan [l4j and Yogindra [15] have extended 
FEM models for anode shape prediction using bit type of 
tools. Ravi Raju [16] has modified model 3GFET~22 for 
cathode design of bare and bare bit tools. In all these 
models however, the stray current zone and the stagnation 
rone have been neglected. 

1 .7.3 BOUNDARY ELEMENT METHOD ( BEM) 

Both the methods, FDM and FEM are domain based 
requiring frequent grid regeneration to compensate for 
the movement of work boundary during simulated ECM, Also, 
both FDM and FEM yield solutions corresponding to internal 
as well as boundary nodes, but in most of the cases, internal 
solutions are unnecessary when predicting the workpiece 
boundary and designing the cathode shape. In BEM unlike FDM 
and FEM, boundary alone is discretized into segments. This 
considerably reduces the amount of input data and its prepar; 
ti on time. However, if required the solution to Internal 
points can also be obtained. 

Narayanan, Hindu j a and Noble [17] have applied BEM, f( 
the anode shape prediction in ECM, Linear and quadratic 
elements were used to represent the boundaries. It was 
concluded that with far less computational effort, a 
comparatively better accuracy is obtained using BEM. However 
for a clear understanding of the process solutions at interns 
nodes are also required. 
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1 ,8 ANODE SHAPE PREDICTION IN STRAY CURRENT ZONE AND 
STAGNATION ZONE 

1 ,8«1 STRAY CURRENT ZONE 

The distribution of electric flux lines during ECD is 
such that some portion of the workpiece which is far away 
from the tool also undergoes dissolution, due to current 
lines acting over this zone (Figure 1e5), This effect is 
called as stray current effect and the zone over which the 
stray current acts is called as stray current zone. Stray 
current attack is undesirable not just because of poor 
dimensional accuracy but also because etched and pitted surfa 
are obtained in this zone, as reported by Hoar and Mears [18 
Rather than predicting the anode shape in stray current zone 
for 3 particular type of electrolyte, the trend in literature 
was observed to be of minimizing the stray current attack by 
choosing an appropriate electrolyte or by modifying the 
properties of a given electrolyte* Boden and Evans [19 j 
showed that the effect of stray current attack could be 
reduced to a certain extent by bringing about passivation 
in the stray current zone by addition of small amount of 
another suitable salt to the electrolyte* They suggested 
that the additive could be carbonate in the case of machinin 
of Iron and Nickie, which would form a sparingly soluble 
reaction product with the workpiece and hence bidng about 
passivation by precipitation of insoluble products in the pi 
of stray current zone, NaC 102 was found suitable for machir 



14 


of cobalt, LaBoda and Mcmillan [20] has suggested the 
use of sodium chlorate as electrolyte for machining of 
any steel workpiece including high alloy steels because 
of its low 'throwing power’ v\^iich ensures that the cutting 
is concentrated in areas close to the cathode, as is 
desired, 

Bannard [21 ] has pointed out that the extent of 
stray current could be controlled to a certain extent by 
use of a " thief-anode" (Figure 1,6) which steals away the 
current flux which 5s not normal to the tool face. The 
disadvantage of this method is the difficulty in the 
manufacturing of the tool, 

1 .3.2 STAGNATION ZONE 

Little literature appears to be available till date, 
discussing about the profile of the spike that is formed 
in ECD of blind holes, Larsson and Muzaffaruddin [22] 
assumed radial current density distribution to determine 
the spike hei^t. The purpose of their analysis was to show 
how the shaping performance of electrolytes is related to 
their electrochemical properties. They used the spike heigh 
as an index to demonstrate the effect of solution properties 
on its shaping performance rather than to predict the actual 
spike profile. The use of a simple tubular drill with an 
uninsulated bore constituted a very sensitive test for shapii 
performance, as the height of the spike produced at the bott 
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of a blind hole is shovm to be dependent on both polarisation 
and current efficiency* They showed that the development of 
a current density dependent voltage at the electrode- 
electrolyte interface, and efficiency of an electrolyte are 
far from negligible in determining the shaping performance 
of an electrolyte. This is particularly so vihen nBchlning is 
in generating mode where there is a wide range of current 
density. They concluded that polarization plays an important 
role in reducing the shaping performance of an electrolyte an< 
should not be ignored in any method used to predict work- 
shapes*, Here, the assumptions made are far reaching, since 
the current distribution at the top of the spike will 
certainly not be radial. Also, the change in electrolyte 
conductivity due to heating effects and the development of 
gas bubbles in the electrolyte has not been taken into accouni 

1 .9 PRESENT WORK 

From the above literature survey, it appears that no 
attempt has been made to determine the effect of stray curreni 
attack on predicted anode shape. Also, no analysis appears 
to be available to find out the profile of the spike that is 
formed in the stagnation zone djring electrochemical drilling 
of blind holes. For a c onprehensive and meaningful analysis c 
anode shape prediction during ECD, the analysis of stray curre 
zone and stagnation zone is very important. The work embodiec 
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in the first part of this thesis has been aimed towards the 
development of anode shape prediction model to analyse the 
effect of stray current attack on the anode profile. The 
anode shape prediction model (93FET-22) developed by 
Yogindra [lb] has been modified by incorporating the stray 
current zone. The model has been named as STZFET-22, This 
model has been extended to determine analytically, the 
profile of spike generated in stagnation zone during ECD, The 
analytical anode profile obtained by STZFET-22 has been compan 
with experimental profile obtained by Jain [3] „ 

The second part of the work has been aimed towards the 
development of a conputer aided method for the tool design ' 
for ECM. The two dimensional cathode (tool) shape prediction 
model for ECD proposed by Ravi Raju [l6] has been modified to 
account for stray current zone and stagnation zone. The anode 
shape prediction model STZFET-22 has been used in this design 
package for electrochemical drilling. Tools have been designe 
for producing the specified work shapes, while machining 
under prescribed machining conditions. Designed tool shapes 
have been compared with experimental tool shapes used by 
Jain [ 3 ] for electrochemical drilling. Design of electrolyte 
supply hole in the cathode has been aimed to minimize the spik 
dimensions, especially the height. This is a case similar to 
the tool design for external shaping by ECM process. The tool 
design model has been named as STZDES-22[i 
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CHAPTER II 

THEORETICAL ANALYSIS 


2.1 INTRODUCTION 

The ability to predict 'the variations in process 
parameters for any given operating conditions is a pre- 
requisite for proper design of tools in ECM, In this 
chapter, the expressions for various process parameters 
such as electrolyte flow velocity, electrolyte conductivity 
current density, temperature distribution etc, have been 
derived for electrochemical drilling. In the present work, 

ECD with outward mode of electrolyte flow has been analysed. 
However, the same analysis can be applied for the case of 
inward mode of electrolyte flow with small changes in the 
computational scheme. To make a comprehensive analysis, all 
the five zones (Figure 2,1 ) namely, stagnation, front, transi- 
tion, side and stray current have been analysed individually 
with suitable assumptions mentioned therein. The stagnation 
zone has been divided into three zones as shown in Figure 
(2,1), The three zones are named as stagnation transition 
(S.T, zone), stagnation side zone (S,S, zone) and stagnation 
stray current zone (S.S.C, zone). 
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2.2 ANALYSIS OF E. ECTROCHEMICAL DRILL IN3 

2.2.1 ELECTROLYTE FLOW VELOCITY (U) 

In the present work, the mode of electrolyte flow 
is outward and separate expressions have to be derived for 
stagnation, front, transition, side and stray current zones. 
Th e a ssumptions made for calculation of electrolyte flow 
velocity are: 

C i) Electrolyte is an ideal fluid, 

(ii) The volume of electrolyte flowing at any section 
per unit time is constant. 

( iii) Cavitation and starvation of the electrolyte are 
neglected. 

( iv) Electrolyte velocity across the flow direction is 
assumed to be uniform 

2 . 2 . 1 . 1 STAG NAT ION ZONE 

As the inner tool corner radius is very small 
electrolyte flow in S.T. zone, S,S, zone and S.S.C. zone 
(Figure 2.1) can be considered to be in axial direction. 

The following analysis to determine electrolyte flow velocity 
holds good in all the three zones, i.e. S.T, zone, S.S. zone 
and S.S.C. zone. Referring to Figure (2,2), 

Ar = (itr^ -n(r^ - y)^) 
or, 

A^ = yC^r^ Y “ Y^) (2.1) 


and 
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U = ^ ^5 (2.2) 

•rt ( 2r^ y - y^) 

where, Q Is discharge of electrolyte/unit time. , 

2.2.1 .2 FRONT ZONE 

In this zone, the electrolyte flow is assumed to be 
in radial direction. Electrolyte flow velocity at any 
radius (Figure 2.3a) from the axis of the tool is given by, 

U = ^ = 2%'? 7 

r 

2.2.1 .3 TRANSITION AND SIDE ZONES 

For calculation of electrolyte flow velocity in 
transition zone, following assumption has been made. If 
the X coordinate of the node at iMaich flow velocity is being 
calculated is less than the outer radius of the tool ( r 2 ) 
then that particular node is considered to be in front zone, 
othervd.se it is considered to be in side zone. In side zone, 
the electrolyte flow is assumed to be in axial direction. 


From Figure (2,3c), 

Aj, = -rt (rg + Y)2 _ 3-2 

or 

= It ( + 2r2 Y) 


and 
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U = « 2 (2„4) 

11 ( + 2r2Y) 

2,2,1 .4 STRAY CLERENT ZOME 

For finding out the velocity of electrolyte, Laplace 
equation has been solved for velocity potential. Here, 
the following assumption has been made. In stray current 
zone at edge AD (Figure 2,4), velocity is perpendicular to AD, 
which means vjfi = 0, i.e.J^ = constant along edge AD, For 

convenience, ^ has been taken to be equal to zero along edge 
AD. 

T 

For a two dimensional case the Laplace equation can 
be writt en as, 



Using simple triangular elements (Figure 2.4), 
the velocity potential (0) has been assumed to vary linearly 
within the element throughout the solution domain. This 
velocity potential (0) can be represented by a polynomial 

0® = a + bx + cy , (2,5) 

The equation ^,5) can also be written in terms of the 
interoolati ng functions, at each node. 
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i.e. J?) = 


2 2 


3"^ 3 


L MJ 


(e) 


( 2 , 6 ) 


Substituting equation (2,6) in equation (2,5), 
the residue becomes 


R® = 


32 ^® 

— ^ 

3x 


■ 2^ , 

9 y 


ByGalerkin’s method. 


// N^R dx dy = 0, i = 1,2,3 


i.e 


-2 2^e 
ff 

a X 


+ 


— dx dy 

ay 


0 


Integrating by parts w,r,t, x for the first term and 


w.r.t, y for the second 

term. 

gives 

3 N, 

X 

+ 

9n. 

. 1 

y 

) dx dy = 

j 

(lf^ 

n 

X 


ly- 

where, and 

are 

the 

direction cosines along 


(2.7) 


axis respectively, on the boundary of the domain. 

Using equation (2,6) and writing the equation (2,7) in 


nBtrix form gives 
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3 N. 

/ 

^ 3 x 





dx dy 


'» normal velocity, 


^ia \ 


dB 


= {V 

n 


( 2 . 8 ) 

T 


1 

' » 



Vn’ 


(O 


'"iB 


V. 


n 


dB. 


For an element, the matrix form of finite element 
bion becomes. 


e. 


5i 

CM 

53 

^21 

^22 

^23 

K 3 I 

^32 

^33 


= // 

^ ’ 3 X 


/ . A 




^2 


^3 


c J 


3 Nj 

ilT 



r 



A 





V. 


n. 


aN. 

3^)dx dy 


( 2 . 9 ) 


( 2 . 10 ) 
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1,2,3. 


1 

The interpolating functions for an element are 


"2 = ~e - ^1^3) + >'(''3 - ^ '^>'1 - =<3) 

and 

"*3 = 275 - ^ 2 ''l) X(Yi-Y2) + Y(X2 - X,) 

Where, 

A® = 5 - Vl) (V 3 -W (^3'^1-Xl''3) • 

( 2 . 11 ) 

After finding out element stiffness matrices [ K]^ 
and normal velocity vector for each element, they are assembled 
for the entire domain to get a global stiffness matrix and a 
global velocity vector respectively. On applying boundary 
conditions, given by equation (2,12), 


3^=0 along solid surfaces 

0 =0 along open edge (2,12) 


a n 


= ~ Velocity in the gap at 

the end of side zone 
( Figure 2,4) . 


the system will yield a set of linear equations, vihich could 
be solved f or 0 by any standard technique. Here, Gauss 
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elimination technique has been used. 

Once the velocity potential (0) is obtained at each 
node, the velocities in stray current zone are calculated 
using equation (2.13). 

U = + v^ ( 2 . 13 ) 

where, 

2 A ® 

- X3) + - xp + gg (X^ - Xg) 

Y ST li I > ■■*»—■ iiwii 1.- I tiiPti w imt rn tmi mmh m . n i*— '*■■■ n i ■>— imum 

2A ® 

0 

and ^ is given by equation (2,11 ) , 

2.2,2 CURRENT DENSITY (J) AND CONDUCTIVITY (K) 

The current density can be determined either from 
Ohm's law or from potential distribution obtained by sol\n*ng 
Laplace equation. 

The application of a particular case depends upon 
the nature of Interelectrode gap ( lEG) and electric potential 
lines. Ohm’s law can be applied only when the current lines 
are straight and normal to the electrode surface. But this is 
not true when electrodes with small radii or complex shapes 
are used. In such cases, the current density calculated from 
potential distribution obtained by solving Laplce’s equation 
gives more accurate results. The current density at any point 
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on the workpiece is calculated as 

J •= K (2,14) 

The electrolyte conductivity, K, is a function of 
electrolyte temperature (T) as well as void fraction 
T5pton [7 ] has proposed equation (2*15) for calculation 
of electrolyte conductivity, which takes into account 
the effect of electrolyte temperature and hydrogen gas 
liberation*. 

K = K. (1 + a * at) (1 - a^,)^ (2,15) 


wh ere, 

t 

is the initial electrolyte conductivity and a 
is the temperature coefficient of electrical conductivity 
of electrolyte* 

In equation (2*15), Hopenfeld and Cole [11J assumed 
a value of void fraction exponent, n = 1 *5 for the case of 
uniform void distribution. For the case of non-uniform 
void distribution, especially where bubbles are concentrated 
near cathode, Thorpe and Zerkle [23] suggested a value of 
n = 2,0, 

In stray current zone and stagnation zone, the effect 
of void fraction on conductivity is neglected, hence in 
these zones, the electrolyte conductivity is a function of 
temperature alone* 
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2.2.3 TEMPERATURE (T) 

Follovdng assumptions have been made in deriving 
the equations for temperature distribution: 

( 1) Total heat generated is only due to Ohmic heating. 

The heat generated by other sources such as friction, 
chemical reactions etc, is negligible, 

( ii) Specific heat and density of the electrolyte remains 
constant • 

(iii) Heat transfer through the electrodes is considered to 
be negligible, i.e, all the heat generated is retained 
by the electrolyte, 

(iv) The flow velocities are such that the frictional 
heating of fluid is negligible,. 


2 .2 o 3 , 1 STAG NAT ION ZONE 

In S.T. zone, S.S. zone and S.S.C, zone (Figure 2,1), 
the electrolyte flovj is assumed to be in axial direction. 

The following analysis for calculation of temperature is appli- 
cable to S.T. zone, S.S. zone and S.S.C. zone. Consider 
an element of length dx at a distance x from the starting of 
S.T. zone (Refer Figure (2*2)),, 


By Ohm’s law, current through this element is given 
by. 


(V -AV) 

y 


♦ 


I 


K 2 - y) dx 


(2.16) 
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Rate of heat generated in the element is given by 


dH = I R 


(V -Av) K 2Tc(r^ -> y) dx 


( 2.17) 


Rate of heat carried away by the electrolyte. 


dH = (5 Pg Cg dT . 

Equating ( 2,1 8) and (2^17) yields 


( 2.1 8 ) 


Q dT = 


(V -Av) K 2Tt(r^ - y)dx 


Using K = Kj ( 1 +aAT), the expression for temper- 
ature distribution becomes, 


dT 


1 ^(T-T.) 


= A’ dx 


( 2 . 1 ^) 


wh ere, 


(V - AY)^ 2.t(r^ - y) 


y Pe Ce 


Integrating equation (2.19) with boundary condition, 
at X " T - T^, gives. 


T - T. 


1 

a 


{[exp a A' (x - - 1 • (2,20) 


2 .2.3.2 FRONT ZONE 

Consider an element of thickness ’dr’ at a distance 
r, from the axis of the tool (Figure 2,3d), 
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Using Ohm's law, the current I through this element , 
Is given by, 


I = 


(V ~AV) K 2r. r dr 
Y 


( 2 . 21 ) 


Rate of heat generated in the element is given by. 


o 

dH 




(V-Ay) K 2Tcr dr 
Y 


(2o22) 


Rate of heat carried away by the electrolyte. 


dH = Q C dT . 
€ 0 


(2.23) 


Equating (2,23) and (2^22) yields 


Q ^ dT = 
e e 


(V K 27rr dr 


Using K = ( 1 +a, AT), the expression for temperature 

distribution becomes, 

1 +aVT' - Tj~) "" ^ ^ (2o24) 

where, 

„ 2TC K. (V - AV)^ 

A = g — 7 ♦ 

The integration of equation (2,24) with boundary 

y 

conditions, at r - ^ ^ gives. 
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T - T: 


= i ^ exp [ A ( ■ 


.)]- 1 > (2,25) 


2.2,3 o3 TRANSITION AND SIDE ZONES 

In these two 7ones, the electrolyte flow is assumed to 
be in axial direction. Consider an element of length dx at 
a distance X from the end of front zone, i,e, from 
(Figure 2.3e), 

By Ohm’ s law. 


I = 


(V - AV) K aitCr^ + Y) dx 

? 


(2.26) 


From equation (2,26), the rate of heat generated in 
the element is given by 


o 

dH 


I^ R 


(V - ^\fr K 27 r( r^ + Y) dx 




‘(2,27) 


Equating equations (2,27) and (2,23) gives. 


, (V -AV)^ K + Y)dx 

Q dT = y ^ 


; 2 . 28 ) 


or 


dT 


1 + a\'^“-”Ti) 


It t 

= A dx 


( 2,29) 


where. 
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2-irK.(V - AV)^ (rp + Y) 

A = £ (2,30) 

Q % Y 

Integrating equation (2,29) with boundary condition 
X = >^ref’ ^ "^ref 

T - ^ tB exp [oa’" (X - ] - 1 1 (2.31 ) 

where, 

B = 1 + - Tj) (2.32) 


2. 2. 3. 4 STRAY OJRREMT ZONE ' 

Consider an element of length ’1’ at a distance r^ 
from the axis of the tool (Figure 2,5), 

By Ohm' s law, 

(V -AV)K 2 Ttr’ dr' 

I = j (2,33) 


From equation (2,33) the rate of heat generated in the 
element is given by, 


-r2 ^ _ (V -AV) K 2 7t r’ dr’ 

dH = I Rg - 


( 2,34) 


Rate of heat carried away by the electrolyte is given 


by, 


dH = Q ^ e dT 


(2,35) 


Equating (2*35) and (2,34) gives. 


dT 


2(v -Ay)^ KTIT' dr» 
1 


(2,36) 


After rearranging the terms. 


dT 


1 + a(T - T^) 


A IV , , , 

= A r ' dr ' 


(2,37) 


where. 


IV 


2 (V -^V)‘ 


1 Q P e 


(2.38) 


Integrating equation (2,37) with boundary condition, 
at r’ as R, T = T^_ gives. 


IV 

T - exp [a I ( r’ ^ - R^) ] - 1 } 


(2.39) 


where, 


B, = 1 +a(T„ - T^) , 


(2.40) 


In calculation of t^peratures in Front, Transition, 

Side and stray current zones, initial temperature is taken 

to be equal to temperature at the end of stagnation zones 

say, T , Similarly, initial conductivity of electrolyte is 
s 

taken to be equal to the conductivity at the end of stagnation 
zone, say, K„ , 
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2.2.4 VOID FRACTION (a^) 

As the electrolyte flows downstream within the machin- 
ing gap, due to electrolysis, hydrogen will evolve at cathode* 

This hydrogen evolution, usually reduces the effective conductivity 
of the electrolyte, so that the local anodic dissolution rate 
varies dovmstream until an equilibrium is achieved. 

Considering a control volume in IBG (Figure 2,6), using 
the principle of conservation of mass, Thorpe and Zerkle [23] 
arrived at the governing differential equation, and the 
expression for void fraction was obtained by solving it as 
f ollows; 


a 


V 


« * 

B X 


1 + B X 


(2,41 ) 


wh ere, 


* 

B = 


crPg 


4, Y 

^ ~ LS rectilinear flow, 


* 

X ^ 


X(r. + X/2) 

— . • >i.^ — for radial flow 

LS r^ 


and 


S = 


V V 

Tipyr- 


Here, E is the electrochemical equivalent and the suffixes > 
a,e and g stand for anode, electrolyte and gas respectively. 


I 


1 
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In the present case, the gas slip rati-o cr in the variable B 
is assumed to be unity, ice,, velocity of gas is same as that 
of velocity of electrolyte* 

The effect of void fraction on conductivity has however 
been neglected in the stray current zone and stagnation zone 
due to the mathematical complexity of the problem* 

?.c2,5 ELECTROLYTE PRESSURE (P) 

The electrolyte pressure [1] required to maintain 
the rate of electrolyte flow in a particular gap is given by 

P = 4 - (2*42) 

where, is the pressure required to over come inertia force 
and P 2 is the pressure required to overcome viscous force o 

The expressions for P^ and P 2 are different for laminar 
and turbulent flows. The type of flow in the lEG can be 
determined from Reynold’s number (Rj-j)» given by, 

R = -S— ^ (2.43) 

" R e 

where, P q viscosity of the electrolyte and D is the 

hydraulic mean diameter v(hich is given by, 

n - 4 X Area of cross section 

” Wetted pierimeter “ 


For a cylindrical tube, D = d, for a rectangular channel of 
depth Y, D = 2Y and for the case in \Ahich fluid is flowing 
through the annulus, D ( d^ d 2 ) . 
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For the laminar flow 2000) with a parabolic 

velocity distribution, the pressures and P 2 are given by 


P 


1 



and 





(2,44) 


for a narrow tube of diameter d. 

The flow is usually turbulent ( > 2000) in ECM, 
From Bernouli’s equation. 




The value of P 2 is given by Blasius equation. 


0,3164 P ^ X 


e e 


2D 


(2,45) 


The pressure P 2 » in turbulent flow is a function of 
Reynold's number. At high Reynold's number inertia force 
will be more than the viscous force which results in 
higher value of P^ « 


2,2.6 FEED RATE ( f ) 

In el ectrochend-cal drilling, as the machining progresses, 
the side zone becomes tapered and it no longer remains a zero 
feed rate case. For the calculation of IB3 in the side zone, 
the feed r^te vector normal to the work surface should be 
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considered, as given by equation (2^46) 

- foCosO (2.46) 

where, 0 is the angle of inclination between the feed direction 
and normal to the work surfa ce , The values of ’9’ for each 
element in the side zone and transition zone are calculated, 
separatel y. 

If and are the interelectrode gaps in side 
zone at points 1 and 2 respectively (Figure 2,7), then the 
angle 9 at the point 2 is given by equation (2.47), 

0 = tan ^ ir? 7^) •* (2,47) 

In the transition zone, the angle Q at any point is 
given by the intersection angle between the normal drawn at 
a point and the tool feed direction (Figure 2,7). 

In stagnation zone analysis, equations (2.-^6) and (2,47) 
have been used for calculating the feed rate in stagnation-side 
zone (Figure 2,1 ), Since stagr^ti on-transition zone is similar 
to transition zone, feed rate in stagnation-transition zone has 
been calculated on the same lines as feed rate in transition 
zone. Unlike the other zones, in stray current zone feed rate 
does not play a role in determining the anode profile and is 
hence ignored. 
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2,2*7 IWERELECTRODE GAP (183) 

Electrochemical reactions would attain an equilibrium 
only after 't' approaches infinity and IBS can be computed 
from equations (2,4 8) and (2,49) while machining with zero 
feed rate and finite feed rate [ 1 ] respectively.; 

Y = (Y^ + 2Ct)^/^ (2.48) 


where, 


C - 


AK ( V - AV) 

Z p F 
w 


and 


t 



Y) 


+ In (^” 




(2,49) 


Equation (2,49) is an implicit equation and can be solved 
by iteration process. Equation (2;i48) predicts an infinite 
value of the gap as t approaches infinity. However, in 
practice, as the lEG increases, the current decreases, so 
that the MRR gradually decreases and for a very large gap 
the process would come to a stand still,; 

For a very small interval of time At, Jain and Pandey [13] 
derived an equation (2,50) for determining lEG which holds good 
for both, zero and finite feed rate . 


= Y +(C' -f)At 
o 


whe re. 


Y 


(2,50) 
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A 



Equation (2^50) is based on the assumption that the 
current density remains constant over a small element of 
length, dX and for a small interval of time, A to This 
equation yields values of Y which are very close to those 
obtained from equations ( 2o48) and (.2o49) and agree well with 
experimental data [3] « In the present work equation (2,50) 
has been used. 


2,2.8 CURREKT EFFICIENCY (h) 

In deriving the expression for metal removal rate by 
Faraday’s laws of electrolysis, current efficiency was assumed 
to be 1005^, But in actual machining process, the MRR depends 
on the current efficiency achieved [ 1 ] , i,e., the proportion 
of total current that is used for the removal of metal from , 
the anode. The current efficiency is given by. 


Actual tnetal removed 
theoretical metal’ 
removed 


X 100% 


(2,51) 


In actual practice, the values of current efficiency 
may be higher or lower than 100% depending on the conditions 
of machining. Current efficiency values less than 100% are 
obtained because of the side reactions at the anode, e,g„, 
evolution of oxygen gas or Fe ions instead of Fe ions 
etc , 
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If the dissolution of the anode occurs at a valency 
lov/er than the one used for calculation, the current 
efficiency obtained may be higher than 100/^o Moreover, 
selective dissolution may occur at the material irregulari- 
ties such as at grain boundaries [24] and hence may cause 
removal of chunks of material o In this case, the actual 
material removed is more than the theoretical one, which 
results in current efficiency higher than 100?^* 





Tr. Zone — Transition Zone 

S.T. Zone —Stagnation - Transition Zon€ 

5.5. C Zone-Stagnation - Stray Current Z 

5.5. Zone — Stagnation - S ide Zone 


Fig. 2.1 Various zones in ECO 




(a) Electrolyte Flow Area 


(b) Current Flow Area 


Fig. 2.2 Electrolyte and current flow area in stagnation 
zone for outward mode of electrolyte flow. 



(b) Front zone (section at A-A) (c) Transition < 

Electrolyte Flow Area 



(d) Front Zone (e) Transition c 

Current Flow Area 


Fig. 2.3 Electrolyte and current flow area in 
transition and side zones. 
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Boundary Conditi'on 
AD 0 = 0 

AB, CD C3O/0n) = 0 

BC {30/00) =Vn 


Fig, 2.4 Finite element discretization & velocity boundary 
conditions in stray current zone. 




Fg 2 5 Current flow area in stray current zone for out 
ward mode of e ectro yte f ow 
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CHAPTER III 
TOOL DESIGN 


3.1 INTRODUCTION 

Electr cxhemical drilling is a process of making holes 
of any shape and size using the principle of controlled 
anodic dissolution. The electrolytic cell for an ECD process 
is illustrated in Figure (3«1), The cathode is simply a 
straight tube, shaped for obtaining the desired work geometry. 
The electrolyte is carried to the drilled hole in the work- 
piece ( anode) through the tool, and metal removal takes place 
as the electrode is fed in the downward direction, 

3.2 PREDICTION OF ANODE SHAPE 

For the purpose of analysis of the process, the entire 
nHChining region has been divided into five zones (.Figure 3,2), 
namely, stagnation, front, transition, side and stray current 
zone. To reduce the complexity in the present problem, the 
stagnation zone has been analysed separately from the other 
four zones. Further, depending upon the location of inlet and 
outlet of electrolyte, the electrolyte flow mode can be class!- 
) fied as outward, inward or side flow. 

In the present work, the finite element technique has 
been applied to the circular hole drilling with bare tool for 
A Qijtward mode of electrolyte flow. Anode shape prediction 

model developed by Jain and Yogindra [15] has been modified to 

■ 



48 


predict the anode shape in all the five zones mentioned 

above. In the anode shape prediction model developed 

here (STZFHT-22), the field distribution is obtained by 

2 2 

solving the Laplace equation = O) with 

ax ay 

the following boundary conditions: 


( 3 ) = 

(ii) 0 = 


0 along cathode (tool) 

(V - AV) along anode (work) 


0 (along open edges). 


Since ECD is a moving boundary problem, the domain 
shape and size will change after each computational machining 
cycle time. At sec. As the tool penetrates in the work surface, 
the total number of nodes would also increase. The coordinates 
of each node, therefore should be recomputed after each compu- 
tational cycle. Figure (3,3) shows the finite element discre- 
tization and renumbering of nodes after one computational cycle 
in all the zones except stagnation zone. Figure (3,4) 
illustrates the same for stagreition zone analysis. Here, an 
assumption has been made that the feed to the tool is given 
only after machining has taken place for a time period of A t 
seconds, i.e. the process is considered as intermittent cutting 
rather than continuous one, as in actual ECM, 

Due to complex nature of ECM process, majority of 
models proposed for anode shape prediction in ECD, do not 
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account for the simultaneous variations in the process 
parameters such as temperature, conduc t.i vity etc. It is 
felt that for an accurate prediction of anode profile, 
the variations in all these process parameters with machining 
time and their effect on anode profile should be consideredo 
These aspects have been included in the present model, 

STZFET-22, Some problems have been attempted using STZFET-22 
and the predicted anode profiles have been compared with 
experimental anode profiles. 

3.3 TOOL (CATHODE) DESIGN FOR ECD 

The problem of designing a tool to produce a required 
work shape is probably the single most important problem in ECM, 
and in practise, methods of design are still empirical even 
to the point of ’trial and error'. The lEG between tool and 
work is a function of various machining parameters. Hence, 
the tool shape is not congruent with the vjork shape in 
the strict sense. The ECM tool design depends essentially 
on the accuracy of the corresponding anode shape prediction 

I 

model. Therefore, for accurate tool-design is ECM, the anode 
shape prediction model must account for the variations in 
process parameters while arriving at an anode shape, 

Reddy [25] has suggested a 'Correction Factor Method' 
for tool design in ECM, In this chapter, a two dimensional 
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model for predicting a cathode configuration which can 
produce specified anode shape under prescribed machining 
conditions has been described. This model has been sugges- 
ted by Ravi Raju [16 j and has been modified here to include 
stray current zone and stagnation zone. The anode shape 
prediction model (STZFEr-22) described earlier has been used 
for tool design in ECD, This model has been suitably 
modified for the prediction of work-shape with a tapered tool. 


3,4 GENERAL DESIGN PROCEDURE 

Tool design procedure for producing the required work 
shape while machining under given conditions is summarized 
below: 

STEP 1 - INITIAL TOOL SHAPE 

The first step is to assume a tool shape depending on 

the required work shape. For majority of the cases, the tool 
shape can be assumed as congruent to the desired work shape. 

STEP 2 “ WORK SHAPE PREDICTION 

The second step is to predict the work shape 
obtainable from the tool shape assumed in the first step. 

Out of the rany models available for anode shape predlotion 
in ECM, the models based on numerical methods like FM and 
FEM, can predict the work profiles reasonably accurately. 

STEP 3 ” error and CORRBOTION 

The third step is to compare the predicted work 

shape mth the required work shape and to obtain deviations 



between them. If the deviation is more than the specified 
tolerance, correction required to minimise this error has to 
be calculated. 

STEP 4 - MODIFICATION 

The fourth step is to modify the present tool shape 
by applying the correction appropriately. Using this modified 
tool, again the anode shape is predicted. Steps 2,3 and 4 are 
repeated till the predicted work shape with the desired tolerance 
is obtained. 

This is an iterative process and the flow chart for 
this process is shown in Figure (3,6), Tool shape is refined 
progressively in each design cycle. Finally, it will give 

1 

the tool shape which produces the desired work shape. 

From the above mentioned procedure, it is evident that 
any inaccuracy associated with the anode shape prediction 
model will be reflected on the designed tool shape • 

3,5 TOOL DESIGN FOR ECD 

Here, an attempt has been made to design bare tool 
for electrochemical drilling. Here, a case of only cylindrical 
hole (straight sided and tapered both) drilling has been 
considered. However, the model can be applied with slight 
changes, to drilling of non-cylindrical holes as well. Using 
above mentioned procedure the tool design model developed by 
Ravi Raju [1 6] has been modified to incorporate the tool design 
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for stray current zone and stagnation zone. Details of the 
above mentioned steps for the present model are as given 
below. 


3.5.1 INITIAL TOOL SHAPE 

Initial tool shape has been assumed as a cylindrical 
one, with its radius equal to 1 mm lesser than the required 
hole radius for the case of tool designed for uniform diameter 
hole drilling. In transition zone, tool corner radius is aghin 
assumed as 1 mm lesser than the required work surface radius 
in the transition zone. This 1 mm value has been chosen 
arbitrarily, and any other value could as well be chosen, 

3.5.2 PREDICTION OF ANODE SHAPE 

The anode shape prediction model discussed in the 
earlier section has been used for computing the anode profile. 
The IE3’s at the beginning of a computational cycle get 
modified in the following cycle. Also, the tool occupies a new 
position after each computational machining cycle. The 
modified coordinates are assigned to the corresponding points 
on tool and workpiece, to obtain the correct work shape at 
the end of a computational cycle,* 

3.5.3 CALCULATION OF ERROR AND CORRECTION 

The coordinate system for the five zones namely 
stagnation, front, transition, side and stray ourrent zones 
is shewn In Figure (3.5). The system of coordinates has been 
SO chosen that the same error calculation principle holds good 
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for all the five zones. The error and correction at each 
node are calculated as followst 

ERROR Cl) = OWS (I) - RWS (I) (3„1) 

CORE (I) = -ERROR (I) (3.2) 

where, 

CORR (I) a: Correction to be made at node. 

OWS (I) = Obtained work shape at node. 

The obtained work shape at any node in front, transition, 
side zones and, S.T, zone and S.S. zone of stagnation zone 
is represented by Y co-ordinate of that node. However, in 
stray current zone and S.S.C, zone of stagnation zone, 

the obtained work shape at any node is given by the X coordinate 
of that node, 

RWS (I) = Required work shape at I node. 

For representing the required work shape at any node 
the same argument given above for obtained work shape in 
different zones holds good. 

The correction need not always be equal to error in 
all design cycles. It can be calculated from the trend of 
the previous design cycle. It may sometimes consume less CPU 
time [l6] . From the second design cycle onwards, the correction 
may be calculated as given below.^ 



^4 


CORRp.tI} 


CORR^_^(I) X ERRORj^ ( I) 
ERRORj^_^(I) - ERRORj^ (I) 


(3„3) 


where, suffix K stands for the present design cycle and 
suffix K~1 stands for the previous design cycleo 

However, for the sake of simplicity, in the present 
work, equat5 on (3,2) has been used. 


3.5.4 MODIFICATION 

The tool that has been used in the previous design 
cycle is modified by applying correction to Y coordinates of 
the tool at different nodes. The new tool shape is given by 
equation (3*4) , , 

NTS (I) = Y(I) + CORR(I) (3,4) 

wh ere » 

^ 5 ^ I) = Y coordinate of new tool shape at the I 
node 

Yi[ I) = Y coordinate of tool in previous cycle. 

Correction at any node ( CORR (I)) may be positive or 
negative depending on the deviation between the obtained work 
shape and required work shape at a particular node, 

Unlike plane parallel electrodes ECM, in ECD a node on 

d to different nodes on tool for different 


the workpiece is expose 
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times. Thus, if a tool shape at a node is changed, it 
will affect the work shape at many other nodes on Its 
way. This in addition to the complex nature of ECM owing to 
electrochemical, thermal and hydrodynamic factors, makes the 
tool designing more difficult in ECD with bare tool. The 
former problem, to some extent, is eliminated in case of 
electrochemical bit drilling* 


I 
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Schematic diagram of ECD with outward 
of electrolyte flow. 
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Fig. 3.5 Error representation in ECD. 
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CHAPTER IV 

results and discussions 

To ascertain the validity of any analytical model, 
it IS essential to compare the analytically obtained results 
with the experimental data. With this basis in mind, in 
this chapter, a ccmparison between the analytical and experi- 
mental anode profiles obtained during ECD with bare tools ha: 
been made. Also, the tools designed for producing the exper 
mentally obtained anode profiles and the tools used durir^ 
experimentation have been compared. However, since experi- 
mental work profile in the stagnation zone is not available 
it could not be compared with the analytically obtained work 
profile in the stagnation zone. Further, because of this, 
part of the tool designed for obtaining a desired work shape 
in stagnation zone could not be ccmpared with tool used duri 
experiraentati on, 

Jain [3] conducted experiments using brass as materia 
for tools and low alloy steel castings and low alloy steel 
forgings as materials for anodes (Table - 1), The machining 
conditions for ECD experiments are given in Table - 2. The 
shape and size of anode profiles were measured by preparing ■ 
cast replica of a drilled hole, and projecting it on a tool 
maker’s microscope. The method by which anode corner radius 
was measured, has been described in detail in reference [3], 
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For analytical computation of anode profile, modified 
equation (2.50) suggested by Jain [3] has been used. 

4.1 AMODE SHAPE PREDICTION 

Using the model SlZFEr~22 discussed earlier, anode 
profiles for different machining conditions have been 
obtained, assuming constant valency of electrochemical 
dissolution of the anode and taking a mean value for the 
efficiency of machining (Table —3). Voltage fluctuation 
history (Table - 4) has also been used while computing 
the anode profiles. This would obviously predict more 
accurate anode profile unlike the case in which mean voltage 
is used, since during experimentation voltage was varying 
rather than renaining constant. The anode profile obtained 
for stagnation zone by model STZFET-22 has been superimposed 
on that obtained for other four zones so as to get the profi 
for all the five zones (Figure 3,2), These profiles for 
different machining conditions have been shown in Figure 
(4.1). 

It can be seen from Figure (4.1) that in majority of 
the cases, the electrochemically drilled holes have the shaf 
of an elongated S(/) , The large overcut at and near the t< 
surface of the workpiece is attributed partly to the stray 
current attack and partly to the fact that this section of i 



63 


workpiece has been subjected to electrochemical dissolution 
throughout the period of ECD, A good quantitative agree- 
ment between theoretically predicted and experimentally 
obtained anode profiles was observed (Figure 4,1 ), The 
small deviations between the analytical and experimental 
anode profiles in front, transition, side and stray current 
zones are due to the following reasons: 

(i) The valency of dissolution plays an important role 
in determining the shape and size of the computed 
anode profiles. The present analytical model assumes 
some constant value of valency of dissolution. But in 
actual machining, the valency of dissolution may be 
greater than the assumed one [26] and hence the 
assumed lower valency would result in larger 
computed overcut. In the present work, valency of 
electrochemical dissolution has been assumed to be 2, 

But as reported in [26] for the conditions prevailing 
in present case of ECD, iron may dissolve with a 
valency of 3 instead of 2. This would reduce the 
overcut and thereby narrow the deviation between the 
experimental and analytical anode profiles. 

The electrolyte conckictivity within the lEG is a 
function of both temperature and void fraction as given 
by equation (2.15), The value of exponent ’n’, in 
equation (2,15), depends on the nature of distribution of 
voids in the For uniform void distribution. 
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Hopenfeld and Cols [11] suggested a value of n = 1 *5,^ 
and for non-uniform void distribution, especially 
when bubbles are concentrated near the cathode, 

Thorpe and Zerkle [23] suggested a value of n = 2,0, 

But in actual machining, it is extremely difficult 
to know the exact void distribution. Owing to this 
reason, an average value of n = 1,75 has been assumed 
for computing void fraction. This approxinBtion 
might also have contributed to some error in calcula- 
tion of overcute. In stray current zone, the void 
fraction has been ignored, to keep the problem 
mathematically tractable, and hence the conductivity 
is assumed to be a function of temperature alone. 

This may also lead to some vari.ation between analytical 
and experimental anode profiles in stray current zone, 

(iii) In some of the tests, due to sparking or some other 
unavoidable reasons, the machining had to be stopped 
and the IBB readjusted. This would change both 
the current and voltage variation history and hence would 
lead to deviations between predicted and experimental 
anode profiles. 

( iv ) In the present analytical model, no consideration 
has been given to the effects of grain size, grain 
boundary attack, selective dissolution of various 
alloy components, and other metallurgical aspects. 

This is expected to affect the value of computed 
overcut 0 
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Since no experimental results are available about 
the spike profile in the stagnation zone, the analytically 
predicted anode profile in this particular zone could not 
be compared. However, the analytically obtained spike 
profile seems to be some what larger than the one that could 
have been normally expected. This discrepancy can be 
assigned to the following reasons: 

(a) Reason (i) menti oned above , 

( b) No consideration for the effect of void fraction on 
conductivity. Conductivity has been assumed to be a 
function of temperature alone, 

(c) Use of simple triangular elements in finite element 
discretization of the stagnation zone which may not 
fit exactly with the work or tool boundary, leading 
to an error. 

4.1 .1 COMPARISON BETWEEN THE MODELS STZFEr-22 and SBFET-II 
From Figure (4.1). it is evident that the model 
STZFET-22 gives better agreement with experimental anode 
profiles in front, transition and side zones, as compared 
to those obtained from SBFET-11, Also, the model SBFET-1 1 
is incapable of predicting the anode profile in the stray 
current zone and stagnation zone. The model SBFET-11 :s a 
one-dimensional model and is based on the assumption that 
the lines of electric flux are straight and normal to the 
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electrode surfaces. The current density in the IB3 Is 
calculated from the Ohm's law. But in complex shaped 
workpiece, especially at curved sections, the electric 
flux lines are not straight (Figure 1,5) and hence 
Laplace equation has to be solved to get accurate potential 
distribution in the domain of interest. In the model 
STZFEr-22, current densities at different nodes are calculated 
from potent-ial d^ stribut'' on obtained by solving Laplace's 
equation. That is why, anode profiles predicted by STZFET-22 
are more accurate than the profiles predicted by SBFET-11 , 

The accuracy of the results obtained by finite element 
technique also depends on the type of elements selected for 
discretization of the solution domain. The model STZFET-22 
uses two dimensional triangular elements, 'Ahich fit better 
in the curved boundaries, as compared to one dimensional 
linear elements used in model SBFET-11. This also contributes 
to the better accuracy of model STZFET-22, However, use of 
isoparametric elements would further increase the accuracy 
of computed results. 


4 1 2 EFFECT OF TEMPERATURE AID VOID FRACTIOM ON CONDUCTIVITY 
Figure (4,5a) shows the variation in temperature, and 

void fraction and conductivity along the electrolyte flow 
path distance. The magnitude of increase In temperature Is 
not much. From Figure (4.5a) it can be seen that the effect 


T 
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of void fraction on electrolyte conductivity is more signi- 
ficant than that of tanperature rise. The electrolyte 
conductivity decreases drastically with the electrolyte flow 
distance. For clarity, the variation in temperature and 
conductivity along the electrolyte flow path in stray 
current zone is shown separately in Figure (4o5b')o ^t 
should be noted that in stagnation zone and stray current 
zone, the effect of void fraction on condjctivity is neglected. 
Hence, in these zones the electrolyte conductivity is a 
function of temperature alone, 

4.2 CATHODE DESIGN 

Figure (4,2) shows a comparison between tool shape 
designed by correction method and the tool shape used during 
experimentation for the case of ECD with bare tools, 
under different machining conditions. Here, the required and 
obtained work shapes are coinciding because the design procedure 
is continued till the deviation between the required and 
obtained work shapes lie within an accuracy cf 10 pm. There ^ 
is a good agreement between the designed and experimental 
tool used under the same machining conditions. Kinks, just 
opposite to the unevenness in the experimental profile, are 
observed in the designed tool shape. This is particularly 
observed at the end of transition zone. 

Figure (4,2) also shows the tool designed by Ravi Raju [16] 
using model DES-22. The tools designed by model STZDES-22 and 
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model DES-22 show a close agreement between the two. 

However, the model DES-22 is incapable of predicting 

the tool shape by taking into consideration the stray 
current effect. 

Figure (4.3) shov;s the designed bare tool (using 
model STZDES-22) for producir^ a straight sided hole with 
top surface of the workpiece square with the hole walls 
(i<,e* the top surface of the workpiece is flat and has 
not undergone any damage due to stray current attack). 

The tool shapes designed by model DES--22 have also been 
plotted for conparison (Figure 4,3), Once agadn, the stray 
current zone has been neglected in model DES-22, 

Figure (4.4) shows the designed tool (w^th an uninsula- 
ted bore) for obtaining the desired spike profile. The desired 
spike profiles for different problems have been shown in 
Figure (4,4), It should be noted that the model STZDES-22 
is capable of designing the electrolyte hole profile for any 
desired spike profile. In the present work, the spike profile 

has been chosen such that it is similar in nature to the spike 

1 

profile obtained by experimental tool, but lesser in height. 
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CHAPTER V 

CONCLUSIONS AND SCOPE FOR FUTURE WORK 

5.1 CONCLUSIONS 

The following conclusions have been drawn from the 
present work: 

- The model STZFET-22 is capable of predicting the 

1 

nature of variation of current density, temperature, 
electrolyte condJct3.vity etc, in all the five zones 
during Electrochemical drilling.. 

Designed tool (cathode) shape is not the comple- 
mentary shape of desired workpiece shape reduced by 
some machining gap dimensions o 

There is a good agreement between the designed tool 
and the tool used during experimentation. Thus 
correction method is capable of designing tools 
reasonably accurately for Electrochemixal drilling. 

The accuracy of correction method deoends upon the 
accuracy of anode shape prediction model. Any 
discrepancy is anode shape prediction model will be 
reflected in the designed tool, 

5,2 SCOPE FOR FUTURE WORK 

The work embodied in this thesis can be extended in 

the following directions; 


i 


p 
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( i) An attempt should be made to improve the accuracy of 
the tool design model by improving the anode shape 
prediction model in the follovdng ways: 

(a) More accurate electrolyte flow velocity distri- 
bution should be obtained by solving Navier- 
Stoke's equations, 

(b ) In the present analysis, effect of void fraction 
on conductivity in stagnation zone and stray ' 
current zone has been neglected. It would be of 
interest to know to what degree this assumption caused 
the predictions to deviate from measured results. 

So effect of void fraction on conductivity must be 
analys ed, 

(ii) Extensive data should be collected about the machining 
efficiency, electrochemical equivalent of different 
metals and alloys under different situations, so that 
accurate prediction of anode profiles can be made, 

( iii) Careful experiments should be conducted to find out 
the spike profile in stagnation zone. The final 
comparison between the experimental and calculated 
results must be left to show whether the assumptions 
made in the analytical model leaves the result suffi- 
ciently legitimate for the purpose in hand, 

(iv) For better accuracy, models based on two and three 
dimensional analysis of ECD using isoparametric and 
higher order elements should be evolved „ 
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(v) For the practical applications (for example, 
elliptical, parabolic shape drilling, etc^) 

the tool design models based on three dimensional 
analysis should be evolved, 

(vi) This tool design procedure can be extended to 
external shaping as welhi. 
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APPEmiX - I 


COMPUTATIOMAL DETAILS 

This appendix describes the development of generalised 
computer programs for the computation of anode profiles obtained 
during ECD using a bare tool (STZFET-22) and for computation 
of designed tools during ECD (STZDE$-22)» In these models, 
three noded triangular elements have been used for discre- 
tizing the solution domaino The potential distribution has 
been obtained by solving the Laplace equation. These programs 
are capable of automatically generating the global stiffness 
iTBtrix, applying the prescribed boundary conditions and storing 
the stiffness matrix in the banded form. The resultant set 
of equations have been solved by Gauss Elimination Technique, 

T 

The computational schemes used in anode shape prediction model 
STZFEr-22 and tool design model STZDES-22 are shown in 
Figures A,1 and A, 2 respectively, 

BRIEF DESCRIPTION OF SUBROUTINES USED 

( I) SUBROUTINE INDATAt This subroutine reads and prints 
the input data for the problem^ 

(II) SUBROUTINE XYCORD: It generates the node numbers, 

X and Y-coordinates for each node, connectivity matrix 
for the finite element mesh, the various arrays and 
assigns the boundary conditions for the nodes on the 


tool and workpiece. 



t III) 


SUBROUTINE BANDW: It determines the band width of the 


( IV) 

(V) 

(VI) 

(VII) 

(VIII) 

(IX) 

(X) 

(XI) 


global stiffness, matrix , 

SUBROUTINE STIMAT: This subroutine generates the 
coefficients of elemental stiffness matrices and 
stores the assembled global stiffness matrix in 
banded form, 

SUBROUTINE ABCOND; This subroutine applies prescribed 
boundary conditions for the nodes on the tool and 
the workpiece, 

SUBROUTINE BANSOL: This solves the set of linear 
simultaneous equations by Gauss elimination method 
and gives the potential distributionc 

SUBROUTINE VELSTI: This subroutine generates the 
elemental stiffness matrices and assemble them is global 
stiffness matrix. It also generates the velocity 
boundary matrix , 

SUBROUTINE ABCVHl.: This subroutine applies the velocity 
boundary conditions, 

SUBROUTINE VELSOL; This subroutine solves the set of 
linear simultaneous equations by Gauss’s elimination 

method, 

SUBROUTINE VELCAL: This subrojtine calculates the 
velocities at various nodes on the workpiece, 

SUBROUTINE TE^iCON: This subroutine calculates the 
temperature and current density at various nodes on 


the workpiece 
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( XII) 

( XIII) 

(XIV) 

(XV) 

(XVI) 

( XVII) 

(XVIII) 

( XIX) 


SUBROUTINE CONfAOD: This subroutine calculates 
the electrolyte pressure and void fraction within 
the IB3, It also modifies the electrolyte conductivity 
by taking into consideration the void fractiono 

SUBROUTINE FHDFIAT: This subroutine determines the feed 
rate at various nodes on the workpiece » 

SUBROUTINE MRR: This subroutine calculates the MRR 
and IBS at different nodes on the workpiece, 

SUBROUTINE MODY; This subroutine generates new elements 
at the end of each computational cycle and calculates 
the Y-coordinates of each node, 

SUBROUTINE MODX: This subroutine generates the X-coor~ 
dinates of each node at the end of each computational 
cycle, 

SUBROUTINE PRTOUT: . This subroutine prints out the 

T 

required results in a specific format, 

SUBROUTINE INITIL: This subroutine initializes the 
global stiffness matrices and the potential distribution 
matrices in each cycle of computation. 

SUBROUTINE DESIGN: This subroutine calculates the 
error and correction and develops the new tool shape. 
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TABLE 1 

Identification Number and Hardness of Work Materials [3 ] 


Work Material 

Ide ntification 
number 

Average 
hardn^s 
( BHN) 

Low alloy steel 

casti ngs 

200-299 

239 

Low alloy steel 

castings 

400-499 

187 

Low alloy steel 

f orgings 

600-699 

21 8 

Low alloy steel 

f orgings 

700-799 

220 


Work Material Composition 


Alloying 

element 

Low alloy 

200-299 

steel castings 

400-499 

C 

0.094 

0.234 

S 

0,018 

0.01 1 

P 

0,020 

0,014 

Si 

0.317 

0.500 

Mn 

0^400 

1 .440 

Ni 

- 

0.095 

Cr 

0.760 

0.260 

Mo 

0.920 

0.073 

Cu 

- 

0.065 

A1 

- 

0.,055 


Va 

Other impurities 


Very small 


5 to 6 
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TABLE 2 

Experimental Conditions used during ECM tests f 3] 


Work mterial 

Tool material 

Tool diameter 

Mode of electrolyte 
s upply 

Applied voltage 

Electrolyte flow 
ra te 

El ectrolyte 
conductivi ty 

Machining time 


M.S, low alloy steel casting and low 
alloy steel forgings 
Brass 

7 mm to 1 8 mm 
Radial outward flow 

! 

10-60 V 

—A -5 3 

6,7x10 " to 15x10 m /sec 

-1 -1 

0e003 to OoOII Ohm mm 
600 to 3600 sec. 
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TABLE 3 


Experimental Conditions used 

durinp ECD 

of Cast 


Low Alloy Steel 

XTl 




f = 0.000571 

mVs ; = 

1 ,0 mm 



Workpiece 

K. 

E 


% 

No. 

( Ohm”* mm”"' ) 

V 

(. volt) 

c 

( ampere) 


T ool“1 

r^ = 4.0 

rtc=1 .245 



413 

0.00725 

17.75 

13.1 0 

90.15 

408 

0,. 007 25 

15.31 

16.24 

91 .70 

409 

0.00725 

21 .84 

14.14 

95.00 

406 

0,00725 

12.94 

1 3.28 

86.90 

Tool-3 

r^ = 4.15 

^c=’ •S'' 



407 

0.0063 

12,93 

15,00 

93.15 

421 

0.0070 

14.89 

16,40 

96.20 

424 

0,,0060 

18.50 

19.86 

92.00 

425 

0.0066 

20,,49 

18.80 

91 .86 

Tool-4 

- 4,505 




429 

0.0060 

11 ,57 

1 8.53 

94,00 

401 

0.00547 

12.975 

22,15 

91 .82 

41 8 

0.00532 

15,95 

18,40 

92.50 

427 

0.00532 

17i62 

- 


402 

0.00600 

15.87 

1 8.40 

93.70 
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TABLE 4 


Voltage Variation during Electro-»cheniical DrillAnq (bare 
type of tool) [3] 


Job No, 

421 

Time in 
min. 


0 

20„5 

1 

20,0 

2 

19.0 

3 

18.0 

4 

18.0 

5 

16.5 

6 

15.5 

7 

1 5,5 

8 

15,0 

9 

15,0 

10 

15.0 

1 1 

14.5 

1 2 

14,0 

1 3 

14,5 

14 

14.5 

15 

14,5 

16 

14,0 

17 

14o0 

18 

13.5 

19 

13.5 

20 

13.5 

21 

13,5 

22 

13.5 

23 

13,5 

24 

13.5 

25 

13,25 

26 

13,0 

27 

1 3,5 

28 

13.25 

29 

13.0 

30 

13,0 


425 

407 

409 

21 .25 

19,;5 

20.75 

21 .0 

1 8,5 

20.25 

20,75 

17.5 

20,0 

20.25 

17.0 

20.0 

20,0 

15.5 

19,8 

20.0 

1 4.5 

1 9 .75 

19,8 

1^,5 

19.8 

19.75 

14.0 

19.8 

19,8 

14.0 

19,6 

19,8 

13,75 

19.25 

19.6 

13.0 

19.0 

19,25 

12.5 

19.0 

19.0 

12,5 

19.0 

19.0 

12.0 

1 8i8 

19,0 

11 .75 

18.5 

18.8 

12.0 

1 8,. 5 

1 8,5 

11 .5 

18..5 

18,5 

11 .5 

18,5 

1 8,5 

11 ,0 

1 8.25 

1 8,,5 

11 ,0 

18,5 

18,5 

11 .0 

1 8.5 

1 8.5 

11 ,0 

1 8,,5 

1 8.5 

1 0.75 

18, .5 

18,5 

10,75 

18, .5 

1 8,5 

11 ,25 

18.5 

1 8,5 

11 ,0 

18,25 

1 8,5 

11 .0 

1 8,25 

18,5 

11,0 

18,25 

1 8,25 

11 .0 

1 8.0 

1 8,25 

11,0 

18.0 

18,0 , 

11 .0 

18,0 


424 

408 

427 

21 ,0 

16,0 

24.0 

20,5 

15.5 

23,0 

19.0 

1 8,0 

22.0 

1 8.0 

17,5 

21 ,0 

1 8.25 

17.5 

20,5 

17.75 

16.5 

19^0 

1 8.0 

16.5 

1 8.75 

1 8.25 

1 6',5 

18,5 

18.5 

16, ,25 

18.5 

1 8,25 

16.0 

17.0 

18.0 

15.75 

17.0 

18,0 

15,40 

16.5 

17.75 

15.0 

1 6,5 

17,5 

15,0 

16.5 

17,5 

15.0 

16.5 

17.5 

14,75 

16.625 

17,0 

14,75 

16.5 

16,875 

15,0 

16.5 

15,875 

15,5 

I6i5 

16.5 

15,5 

16,5 

16.5 

16.0 

1 6, .5 

16.0 

17.0 

16.5 

16,5 

17.0 

1 6 .5 

15.5 

16.0 

16.5 

1 5 .5 

1 6.0 

16.5 

15,5 

13,0 

1 6.0 

15,5 

13,0 

I64O 

15.5 

13,0 

15.75 

15,5 

13.0 

15.75 

15.5 

13.0 

15.5 

15.5 

13.0 

15,5 
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TABLE 5 


Experimental Observations about Overcut in Side Zone durincj 

'T§1 ^ 


Job NOe 

425 

427 

421 

4C9 

408 

406 

Depth 

X,jnn 



Overcut, 

mm. 



0,0 

. 2,26 

2,32 

1.71 

2,45 

2.634 

1 .902 

0,5 

1 ,76 

1 .73 

1 ,54 

2.01 

2.46 

1 .52 • 

1 ,0 

1 .695 

1.65 

1,56 

1 ,34 

2.214 

1.75 

1 .5 

1 .61 

1 ,63 

1 .4 

1.7 

2,14 

1 .33 

2,0 

1,6 

1 ,44 

1 .45 

1 .67 

2.02 

1 .36 

2,5 

1 ,51 

1 .51 

1 .31 

1 .57 

1 .964 

1 ,21 

3.0 

1 .52 

1 ,34 

1 .36 

1 .52 

1 ,85 

1.25 

3.5 

— 

1 .42 

- 1.13 

- 

1.76 

1.11 

o 

« 

1,43 

1 ,12 

1 ,25 

1 .42 

1 .67 

1,154 

4,5 

1 ,33 

1 ,31 

1 ,02 

1 ,38 

1 .57 

1,03 

5.0 

1 ,34 

1.02 

- 

1 .34 

1 ,48 

1 ,05 

5,5 

1.23 

1,18 

0,.o91 

— 

1.34 

0.924 

6,0 

1 .11 

0.91 

0,92 

1.26 

1 .28 

0,96 

6.5 

1 .1 

1 ,03 

0.71 

1 .26 

1 .17 

0,84 

7.0 

0.91 

0,61 

0.715 

1 .16 

1 .134 

0,87 

7,5 

0.84 

0.61 


1 .11 

0.95 

0,74 


